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An ordered hexagonal array of aggregates can form in thin ferrofluid layers when an external 
magnetic field is applied. Using the Helmholtz free energy for this system, we calculate the optimum 
spacing for these aggregates. Results show excellent agreement with experimental findings as a 
function of field strength and layer thickness. Our analysis yields a crossover in the exponent for 
the scaling behavior of the aggregate spacing as a function of plate separation, in agreement with 
experiment. For the first time, we report a similar crossover in the scaling behavior of the aggregate 
spacing as a function of the magnetic field. The mechanisms responsible for both crossovers are 
introduced and discussed. 

PACS numbers: 61.20.Gy, 47.54.+r, 82.70.Dd, 75.50.Mm 



I. INTRODUCTION 



Ferrofluids are composed of small single-domain magnetic particles dissolved in a liquid carrier. With no external 
magnetic field, Brownian motion keeps the particles homogeneously distributed throughout the liquid. When a sample 
of ferrofluid is conflned between two glass plates with a separation L, and an external magnetic field Hq is applied 
normal to the plates, the particles aggregate. These aggregates are of approximately uniform size, and the separation 
of the aggregates is periodic with an average spacing d. A view normal to the plates displays a hexagonal pattern of 
aggregates [llllll. 

In this paper we present a new model for predicting the aggregate spacing, and compare the results of this model to 
experiment. We demonstrate that this model shows good agreement with experimental data for the aggregate spacing 
as a function of field Hq and plate separation L. We compare our results to experimental data from Refs. P, fl. 

The model also yields the scaling behavior of the aggregate spacing d as a function of Hq and L. In 1994, Wang 
et al. 4] showed that there was a crossover in the exponent for c? as a function of L. They attributed this crossover 
to an experimentally observed structural change in the aggregates. In 1995, Liu et al. also showed a crossover in 
this exponent for ferrofluid emulsions, but no explanation for the crossover was provided. In our model, this crossover 
takes place, and is not attributed to the mechanism in Ref. i\. We also report, for the flrst time, a crossover in the 
scaling behavior for d as a function of Hq . The mechanisms responsible for both crossovers in our model are identified 
and discussed. 



II. HELMHOLTZ FREE ENERGY 



Experiments show that, when the field is first turned on, the particles form single-particle-diameter chains; the 
system then evolves by the chains sticking together to form aggregates. The time scale associated with the chain 
formation is much smaller than the time scale observed for the formation of aggregates from chains 0. Due to this 
observation, we treat the formation of the aggregates as a quasiequilibrium process. Other examples of this treatment 
can be seen in Refs. ^|lii6{. With this assumption, the separation of the aggregates can be determined by minimizing 
the Helmholtz free energy for the ferrofluid system. 

There are three contributions to the Helmholtz free energy that must be considered: the magnetic energy, the 
entropy and the surface energy. For the purpose of notational simplicity, we introduce the constants a and /3, such 
that 



d = pb, (3 = 




2 



where A'o is the total number of aggregates in the system. Here, b is the aggregate radius, (f> is the volume fraction of 
the ferrofluid, 7 is the packing fraction of the particles in an aggregate and i'^ is the area of one glass plate. These 
relationships are purely geometric, using the assumption that the aggregates are of uniform size and spacing. 

A. Magnetic Energy 

The magnetic energy contains three parts corresponding to the self-energy of the aggregates, the aggregate-aggregate 
interactions, and the aggregate-external field interactions. The model approximates the aggregates as cylinders that 
are uniformly magnetized. 

The approximation of cylindrical aggregates is consistent with the experiments in Refs. |^, Q , where it was found 
that a few particles away from the end of the aggregate the shape was cylindrical. Since we are dealing with L ^ 10 /im 
and 6 ~ 5 nm, the length of an aggregate is on the order of a thousand particles. Given these experimental features, 
the approximation of cylindrical aggregates is a good one. 

The approximation that the magnetization is uniform throughout the aggregates is reasonable. Near the ends of 
the aggregates one would expect the average magnetic moment of the particles to be different than that of the bulk. 
However, the magnetization of the aggregates is determined mainly by the particles in the bulk. A correction for 
surface effects is included in the surface energy, discussed below. The uniform magnetization approximation has also 
been used in Refs. 4, 5]. 

From magnetostatics we know that the effective surface and volume charge densities are am — M -fi and pm = V • Af 
respectively, where fi is the unit outward normal, and M represents the magnetization of an aggregate. For the case 
of uniform magnetization parallel to the external magnetic field, am = M on the top and bottom of the aggregates, 
am = on the sides and pm ~ everywhere. Thus, the total magnetic energy of the system can be calculated as 
interactions between disks of uniform equivalent charge density am = M. 

For the first contribution to the magnetic energy we consider the energy required to accumulate a single disk of 
equivalent charge. This self-energy of the disk can be calculated exactly in the following simple way. Consider a disk 
of radius r <C 6 and uniform charge density am- The potential at a point on the rim of the disk is 

V{r,am) = 4a-,„r. 

Now bring in a ring of uniform charge density am, radius r and thickness dr. The additional self-energy of the disk 
when the ring is added is 

dEseif, = V{r,am)dq = Sna^^r^dr. 
Integrating from r = to r = 6 yields the total self-energy of one disk, 

Eseifi = -namb^. 

To obtain the self-energy for all of the disks in the system, we must multiply this expression by 2Nq and use the fact 
that am = M. The factor of two is due to the fact that each aggregate is replaced by two disks of equivalent charge, 
so there are two layers of A^o disks to consider. The final expression for the self-energy of the system is then 

Eself = —TTaM -. 

This energy partially captures the particle-particle interactions within an aggregate. An additional contribution to 
the single aggregate energy arises from the interaction between the pair of disks associated with the same aggregate, 
as discussed below. 

For the other contributions to the magnetic energy, the disks can be treated as point charges. This approximation 
is valid since the aggregate spacing and length are large when compared to the aggregate diameter (i.e. d ^ b and 
L ^ b). The condition that d 3> 6 is satisfied by having a low volume fraction (j). (The experiments used for 
comparison with our model are done for (j> < 0.20.) To satisfy the second condition, we use experimental observations 
which indicate that larger values of L produce a smaller 6/L ratio. Once the energetically favorable radii have been 
determined by minimizing the free energy, points where 6/L <C 1 are considered accurate and other data points where 
b/L ^ 1 are discarded. 

Using the point charge approximation, the aggregates behave as dipoles, and the interaction energy of the aggregates 
with the external magnetic field is 



Eiig = Noi-m ■ Ho) = -arrLAIHo, 



(1) 
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where m is the magnetic moment of an aggregate. This term does not depend on d, and therefore will be considered 
a constant in the minimization process. 

The last contribution to the magnetic energy is due to the interactions between the disks. We must account for 
interactions between disks in the same plane (i?samc) and disks in opposite planes. For the interaction of the disks 
on opposite planes, two terms must be included. One is the interaction of disks due to the same aggregate (£^opp; i.e. 
one disk directly over the other disk), and the other is the interaction of disks positioned diagonally to each other 
(i?diag)- Using the point charge approximation, as previously mentioned, the total disk interaction energy will be 

/ (P d 

It should be noted that -Eopp arises from particle-particle interactions in a single aggregate, while £^samo and -Ediag 
arise from the interaction between two separate aggregates. 

In Ref. Hong and collaborators found that the magnetization of their ferrofluid sample was well approximated 
by to the Langevin function. Since this functional form for the magnetization is also in agreement with Refs. |^|^, 
and for small (j) in Ref. 0, we have used it for our model. Combining these contributions yields a total magnetic 
energy of 



where 




6d2 



(3^Vd^+L^ 



a-K^M^ + const, (2) 



M=^Ms(coth77-i ), ri=J^Ho. (3) 

Here, fj, is the magnetic moment of the particles in the ferrofluid and Ms is the saturation magnetization of the fluid. 
In the case of a ferrofluid emulsion, each drop has a magnetization described by the Langevin function. In this case, 
the magnetization of an aggregate is still given by Eq. (|3J), with an important reminder that /i is the magnetic moment 
of the particles in the drop, not the magnetic moment of the drop itself. 

Performing an expansion of Eq. (j^J in powers oid/L (with d — f3b), and keeping up to quadratic terms, this 
magnetic energy reduces to the one used in Refs. However, one should note that the assumption c? ^ L is a 

much stronger assumption than b <^ L. For very small volume fractions, one can have the situation where d ^ L with 
the condition b <^ L easily satisfied, since /3 increases as the volume fraction decreases. 

This magnetic energy offers a significant advantage over the model presented in Ref. "iue to the elimination 
of all free parameters in the magnetic energy. With no free parameters in the magnetic energy, we can compare the 
orders of magnitude of the magnetic energy and entropy, as discussed in the next section. 



B. Entropy and Surface Energy 

To calculate the entropy we treat the aggregates as distinguishable. The number of states accessible to the system 
is NqI, and the configurational entropy of the system is S = fcBln(AfoO- For the experimental systems that we 
are modelling, E^^^g/TS ~ 10^, where T is taken to be room temperature. Thus, the entropy is negligible for this 
particular set of experimental conditions. 

The last contribution to the energy to be considered is the surface energy. We have already accounted for particle- 
particle interactions in the bulk in calculating the magnetic energy. The surface energy allows us to include interactions 
on the surface of the aggregates, such as those between the surfactant and the solvent. 

Assuming a surface tension that is constant over the surface of the aggregate (as in Ref. foj), the surface energy of 
the system can be calculated as 

E'surf — 27ra/3(T ( — | -f const, 

V) 

where a is the surface tension of the sides of the aggregate, and the surface energy due to the ends of the aggregate 
is independent of d. There are several ways of dealing with the surface tension (see for example Refs. [ill ll^ V 
We tried each surface tension and found that our results were closest to experiment using the surface energy found 
experimentally in Ref. 'll]. Here it was determined that the surface tension goes as rf/^. Since the surface tension 
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includes interactions between the particles on the surface of the aggregate, we would expect the surface tension to 
saturate at large field values. To build this saturation into our model we assume that the surface energy goes as M'*/^. 
The final form of the surface energy is then 

E,^,f = 2^a/3S^M4/5, (4) 
a 

where E is a constant determining the magnitude of the surface tension. 



C. Total Helmholtz Free Energy 

Combining these contributions, the Helmholtz free energy of the ferrofluid system is given by 



mag 



£^surf — TS ~ E, 



mag 



E. 



'surf 1 



(5) 



with Eyna,g and Eguri given by Eqs. (O and respectively. Factoring out common terms and writing the magnetization 
in terms of (3 we obtain 



F : 




2(3 A L 
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an'^Mg + const. 



(6) 



where 



M 

Ms 



2tt 



cothr^ 



Aee 



M. 
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We minimize this free energy with respect to d numerically using Mathematica 4.0. The values of (3, L, and a are 
fixed by the experimental conditions, so A is the only free parameter in the model. The value of Ms is used only to 
calculate other quantities such as /i and (j). (3 can also be varied slightly, but only within reasonable limits established 
by the value of 4> which is given for each experiment. 

This free energy can be compared to the theoretical models in Refs. Our new model differs from these 

previous models in three significant ways. First, this model has only one free parameter. All other parameters are 
specified by experimental conditions; even this parameter could be determined by experiment. The second difference 
is our treatment of the surface energy using the results of Ref. . Thirdly, this model permits direct comparison 
between the magnitudes of the magnetic energy and entropy. 



III. RESULTS OF THEORETICAL MODEL AND CROSSOVER MECHANISMS 



For the figures below, all parameter values except A are assigned according to the experimental conditions. Then 
A is chosen so that the theoretical results give the best fit to the experimental data. For all of the sets of data, 
T = 300.0 K and ^ = 3.0 x 10~^^ erg/Oe (~10 nm magnetite particles). For Hong et al. P] the parameter values 
are: (3 = 8.0 (calculated from cj) = 0.04 and 7 = 0.71), and A = 3.09 x 10-"* (cm30e)Verg. For Wang et al. |1| the 
parameter values are: (3 = 4.5 (calculated from (j) = 0.12 and 7 = 0.67), and A = 2.69 x 10"'' (cm'^Oe)^/erg. As a 
point of interest, due to the fact that the surface tension includes the interaction of the particles at the surface of an 
aggregate, one might expect the surface tension of an aggregate in a ferrofluid emulsion to be much smaller than that 
of an aggregate in a ferrofluid. Using our parameters, we obtain a ~ 10^"^ erg/cm^ for the ferrofluid emulsion in Ref. 
m and (T ~ 10 erg/cm^ for the ferrofluid in Refs. confirming this expectation. 

We first compare our model to the experiments in Ref. "Tl , where experimental data was recorded for d as a function 
of L for fixed iJp — 300 Oe. They report an exponent of 0.47, changing to an exponent of 0.67 at L « 20 /im. In Fig. 
^ the results of our model are shown as diamonds, squares and triangles for three different values of _ffo. The solid 
lines are power law fits of our numerical results. Although, we predict lower exponents than the data suggests, the 
crossover in the exponent is present. We have included results for three values of Hq to show the effect of changing 
Hq in our model. It is clear that increasing the field produces three significant effects: decreasing the critical value of 
L where the crossover takes place, increasing the value of both exponents, and decreasing the difference between the 
exponents. It is our hope that these effects will be probed experimentally. 

A closer look at our model reveals the mechanism responsible for the crossover. In Ref. ^, the crossover in the 
exponent was thought to be the result of a structural change in the aggregates for increasing L. In their experiment. 
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as L is increased, the aggregates begin to lose their circular cross sections, forming branching structures between 
aggregates. A theoretical model was developed to account for these branched structures, predicting a larger exponent 
as the branching occurs. We have obtained a crossover in the exponent assuming cylindrically shaped aggregates with 
no structural changes, and thus the crossover in our model cannot be attributed to the mechanism in Ref. 

Examination of Eq. lO for L ^ d provides a mechanism for the crossover. In this limit, the free energy reduces to 
the form 

FL^d~Cod + Ci^, (7) 

where Co and Ci arc constants with respect to d and L. This form of the free energy predicts that d oc \IL. When 
d ^ L the exponent will be less than 0.5, with the exponent depending upon the value of /3. Thus, the mechanism for 
the crossover in our model is a competition between terms in the magnetic energy. For L 3> rf, the magnetic energy 
is dominated by the self-energy of the disks and the energy due to disks in the same plane. As L is reduced, the 
contributions due to the disks in opposite planes become more important, until at some critical value of L, they can 
no longer be ignored, lowering the value of the exponent. 

Fig.lHshows the aggregate spacing as a function of H^. The crosses represent the experimental data for L = 10.0 iim 
obtained from Ref. |l|, and the diamonds, squares and triangles represent the numerical results for three different 
values of L. The solid lines are power law fits to our results. More experimental data points would prove beneficial 
in testing the theory in this case; however the theory agrees very well with this experimental data. Fig. |21 shows that 
our model predicts a crossover in the exponent for d as a function of Hq. This is the first time such a prediction has 
been reported in the literature. Three values of L are shown to demonstrate another prediction; the critical value of 
Hq where the crossover occurs, and the value of the exponents, have essentially no dependence on L. 

From Fig. 13 it is clear that, for large L, the crossover in the exponent for d as a function of Hq occurs with 
approximately the same exponents as for small L. Therefore, to understand the mechanism responsible for this 
crossover, it is beneficial to look again at Eq. © in the limit of i ^ d. In this limit, 

FL»d - C2dM^ + C^^M^'\ (8) 

where C2 and C3 are constants with respect to M and d. This predicts that d oc M~^/^. Thus, for small 77, the 

—3/5 

aggregate spacing goes as , and for large 77 the magnetization saturates and the aggregate spacing does not 

depend upon Hq. For the results presented in Fig. [21 this saturation occurs ^ 2000 Oe. Since the magnetization of an 
aggregate is proportional to the Langevin function, it is clear that a change in the magnetic moment of the particles 
in the ferrofluid changes the value of the field where the crossover takes place. Thus, the mechanism for the crossover 
in the exponent for d as a function of Hq in our model is the field dependence of the magnetization of the aggregates. 
This is a direct result of the form chosen for the field dependence of the surface and magnetic energies in our model, 
so experimental studies of this crossover could be used to determine the appropriateness of these field dependencies. 



IV. CONCLUSION 



We have shown that aggregate spacing can be determined by minimization of the Helmholtz free energy. We have 
developed a simple model and compared its predictions to experiments in Refs. Our model predicts trends 

as a function of external field strength and plate separation. The results of our model are in good agreement with 
experiment. 

The scaling behavior of the aggregate spacing as a function of plate separation and external field is discussed. For 
aggregate spacing as a function of plate separation, there is a predicted crossover in the value of the exponent at 
some critical plate separation. This feature of our model is in agreement with Refs. 0,01, where similar results were 
reported. The mechanism responsible for this crossover in our model is a competition between terms in the magnetic 
energy. 

Our model also predicts a crossover in the exponent as a function of external field. This has not been reported 
in the literature to date. The mechanism responsible for this crossover in our model is the field dependence of the 
magnetization of the aggregates. Different forms of the magnetization, surface energy or magnetic energy would 
change this behavior. This result suggests that experiments measuring this crossover could be used to probe the 
forms of these quantities. 
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FIG. 1: Aggregate spacing as a function of plate separation for Hq = 100.0 Oe (upper), Hq = 300.0 Oe (middle) and Hq = 600.0 
Oe (lower), with other parameters chosen to model the system in Ref. 0|. The exponents obtained by a power law fit of our 
numerical results are shown in the figure. 




FIG. 2: Aggregate separation as a function of external field for L = 100 ^ra (upper), L = 50 /im (middle) and L = 10 /im 
(lower), with other parameters chosen to model the system in Ref. The exponents obtained by a power law fit of our 
numerical results are shown in the figure. 



